Coexisting Non-Equilibrium Condensates with Long-Range Spatial Coherence in 

Semiconductor Microcavities 
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Real and momentum space spectrally resolved images of microcavity polariton emission in the 
regime of condensation are investigated under non resonant excitation using a laser source with re- 
duced intensity fluctuations on the timescale of the exciton lifetime. We observe that the polariton 
emission consists of many macroscopically occupied modes. Lower energy modes are strongly local- 
ized by the photonic potential disorder on a scale of few microns. Higher energy modes have finite 
k- vectors and are delocalized over 10-15 /im. All the modes exhibit long range spatial coherence 
comparable to their size. We provide a theoretical model describing the behavior of the system with 
the results of the simulations in good agreement with the experimental observations. We show that 
the multimode emission of the polariton condensate is a result of its nonequilibrium character, the 
interaction with the local photonic potential and the reduced intensity fluctuations of the excitation 
laser. 



I. INTRODUCTION 



There is significant contemporary interest in the study 
of strongly coupled semiconductor microcavities, where 
mixed exciton-photon quasi-particles, two-dimensional 
(2D) polaritons with very small effective mass can be 
created [H, [3, @]- As a result, high density macroscopi- 
cally occupied polariton states can be achieved at high 
temperatures and relatively small optical excitation den- 
sities. Most notably, the formation of macroscopically oc- 
cupied states arising from polariton Bose-Einstein Con- 
densation (BEC), has been recently observed in CdTe 
H, GaAs H 0, and GaN @, |jj microcavities under 
conditions of non-resonant excitation. Polariton conden- 
sates exhibit the characteristic properties of BEC, such as 
long range spatial coherence [4], suppresion of the pho- 
ton bunching effect at threshold Q , long coherence time 
for both the first (gl) and second (g2) order correlation 
functions [10[ and the quantization of vortices [ll| . Al- 
though there is evidence for equilibrated polariton distri- 
butions 0, Qjl , and macroscopic occupation of the k=0 
ground state, polariton condensates are far from thermo- 
dynamical equilibrium: due to the rather small polariton 
lifetimes these states originate from a dynamical balance 
of pumping and losses in the system. There are several 
approaches based on both microscopic calculations and 
generalization of the Gross-Pitaevskii equation, which de- 
scribe the consequences of the non-equilibrium aspect of 
the interacting polariton BEC system [l4L [TBI, ITR Il7j]. 

The imperfections of semiconductor microcavity struc- 
tures result in a disordered photonic potential landscape, 
which plays an important role in polariton propagation 
and especially, formation of 2D macroscopically occupied 
polariton states. The influence of photonic disorder was 
first discussed for the case of the optical parametric os- 
cillator [18[, where macroscopically occupied signal and 



idler states are formed due to direct scattering of the res- 
onantly excited pump polaritons. It was shown that the 
photonic disorder strongly affects the real space distribu- 
tion of the signal [19], which consists of spatially local- 
ized modes emitting at different energies [20j. The non- 
equilibrium polariton BEC in CdTe microcavities excited 
non-resonantly was also observed to consist of several lo- 
calized maxima [4|, which arise from polariton trapping 
by the potential disorder. In Ref. [2l|, it was demon- 
strated that synchronization between condensates local- 
ized in neighboring minima is possible. In addition, it was 
shown that polariton condensates subject to a disorder 
potential may spontaneously exhibit quantized vortices 
due to interplay between the disorder potential and the 
flow and decay of the condensate [111 ]. 

Recent experiments revealing very long coherence 
times of the polariton BEC using CW diode laser excita- 
tion in CdTe structures showed that microcavity polari- 
ton emission consists of multiple narrow modes strongly 
overlapping in real space [10]. Such an observation be- 
came possible, since the diode laser has reduced inten- 
sity noise on a timescale of 100-1000 ps and results in 
reduced fluctuations of the number of excitons and po- 
laritons, which otherwise lead to marked broadening of 
BEC spectra and prevent the observation of the coexist- 
ing condensates jlQj . 

In the present work, we present spectrally resolved real 
and momentum space images of the multimode polariton 
condensate and probe the properties of the spectrally re- 
solved long range spatial coherence of each mode. We 
observe that above threshold condensation occurs into 
several polariton levels. This is a result of the non- 
equilibrium aspects of the system, which prevents ther- 
modynamical equilibrium being reached for the states 
near the bottom of the lower polariton (LP) branch. The 
low energy condensate states are found to be strongly lo- 
calized within a deep potential minimum, whereas the 



2 



higher energy states are spatially extended over few po- 
tential minima. In momentum space, these delocalised 
states typically consist of a few bright maxima dis- 
tributed on a ring. Similar to the phenomenon of co- 
herent back scattering of resonantly excited disordered 
systems, a peak at some k is often accompanied by a 
weaker satellite at — k. 

After the discussion of the experimental observations 
in Sec. HI we present in Sec. IIIII a theoretical analy- 
sis of a multiple frequency condensate based on a mean 
field model for polariton condensation P3], taking into 
account the finite life time and the disorder acting on 
the polaritons. Very good qualitative agreement be- 
tween theory and experiment is observed. Conclusions 
are drawn in Sec. IIVI 



II. EXPERIMENT 

The sample employed here is the one used in Ref. [4] 
where a BEC with extended spatial coherence was re- 
ported. The sample was cooled to 10°K. To avoid heat- 
ing of the sample, quasi- CW non-resonant excitation was 
employed by means of two CW diode lasers at 685 nm us- 
ing a mechanical optical chopper with a frequency of 300 
Hz. The excitation conditions are similar to those em- 
ployed in Ref [l0[| , where multimode emission of polariton 
BEC was observed. The size of the excitation spot was 
about 20 jam. High resolution sub micron imaging was 
achieved by a high numerical aperture microscope objec- 
tive (N.A=0.5). Spectrally and spatially resolved images 
were recorded using a double monochromator of about 
30 fieV resolution and a CCD camera. Interferometric 
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measurements were performed by means of a modified 
Michelson interferometer in a mirror-retroreflector con- 
figuration with active stabilization 0, [ll|, EH . 



A. Spectrally resolved images in real and 
momentum space 

At low excitation powers below the condensation 
threshold polariton emission is quite broad (FWHM 1 — 
1.5 meV) with no resolved spectral features, as reported 
previously 0, [TO]. By contrast, above condensation 
threshold strong polariton emission arises, which consists 
of multiple narrow lines [lOj- FigOJa),(b) show the im- 
ages (dispersions) of polariton condensate emission above 
the condensation threshold in energy- momentum (E-k x ) 
space taken at momentum k y = for P = 1.4 mW 
« 2P t h excitation power, where P t h is the condensation 
threshold. The spectra are recorded for two perpendic- 
ular linear polarizations X and Y which are respectively 
parallel and perpendicular to the crystallographic axis 
(1, 1,0) of the sample. The polarization degree is of the 
order of 80 percent with the strongest emission in X- 
polarisation, which is consistent with our previous obser- 
vation [4]. FigHtc) shows the corresponding spectrum 
obtained by collecting light from the 20 /im excitation 
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FIG. 1: (a), (b) Dispersions of the lower polariton branch for 
two perpendicular polarizations X and Y.(c)The real space 
intensity of the luminescence averaged over mode size. 



FIG. 2: (a) real space spectrally averaged luminescence image 
for X polarization. (b)-(f)local spectra obtained from small 
regions in real space at positions shown by the arrows in panel 
(a). 
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FIG. 3: (a)-(f): state 1-state 6 momentum space imaging at 
the energy of each mode. 
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FIG. 4: (a)-(f): state 1-state 6 real space imaging at the 
energy of each mode. 



area. We are able to distinguish 6 narrow peaks (states 
1-6) in the polariton emission of each polarization. The 
energies of these peaks are separated by 100 — 250/ieV. 
A splitting of w 25/ieV between the X and Y polarized 
narrow lines is observed for each of the modes. Such 
a splitting also exists for the LP mode below threshold 
and is probably due to intrinsic anisotropy in the micro- 
cavity mirrors, leading to a cavity birefringence. Only 
the first 4 lower energy peaks (states 1-4) are observed 
in FigQJa),(b) in (E-k x ) space, since the higher energy 
peaks (states 5-6) have maxima of intensity at k y > 0. 
However, these peaks are clearly revealed in 2D images 
in momentum space as shown in Figj3j 

In order to reveal further information on the mode 
structure without complication of spatial averaging as 
in FigJTJc), spectra were recorded with spatial selection 
from specific regions of ~ 1/im 2 size across the excitation 
spot. The real space spectrally averaged image of the 
condensate is shown in FigJSfa). Fig|5fb)-(f) shows X- 
polarized condensed polariton emission spectra obtained 
for detection at each of the points labeled A-E in FigjSfa). 
Different peaks dominate depending on the detection 
area across the spot. However, as seen in FigJ2(b)-(f) 
the emission from small regions always consists of about 
3 peaks, which suggests that the modes of the condensate 
strongly overlap in space. 

In both polarizations the peaks are narrow: their 
linewidths (FWHM) lie in the range from 150 jaeV down 
to 50 fieV (resolution limited) for the most intense states 
3 and 4 [12]. Such an observation is consistent with 
our previous studies, where long coherence times up to 
200 ps (FWHM= 10 /ieV) were reported for a single 
mode of microcavity polariton emission above conden- 
sation threshold. Strong spectral narrowing and super- 
linear increase of the intensity of each single mode is ob- 
served with increasing excitation power [lOj], which in- 
dicates the build-up of coexisting macroscopically occu- 
pied states (polariton condensates) and condensation into 
both non-degenerate X and Y polarized states. 

The number of condensates observed varies from 1 up 
to 6 — 7, and this, as well as the energy separation be- 
tween the condensed modes, depends strongly on the 
excitation position across the sample. As we show in 
this manuscript the multimode emission above threshold 
arises from interaction of polaritons with the local pho- 
tonic disorder potential. 

To further reveal the origin of formation of the spec- 
trally narrow coexisting condensates we recorded spec- 
trally and spatially resolved 2D images in real and mo- 
mentum space. The images were obtained above thresh- 
old for each of the narrow peaks in FigOJc). Although 
the Y-polarized modes are spectrally distinct from the 
corresponding X-polarized modes, they are observed to 
have an almost identical pattern in both real and momen- 
tum space. Such a similarity arises from the fact that the 
photonic potential, which determines the new polariton 
condensed states, is given by the spatial fluctuations of 
optical thickness of the cavity due to defects in the dis- 
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FIG. 5: (a)-(f)real space imaging of the interference pattern 
of modes state 1-state 6. 



tributed Bragg reflectors and thus is identical for both 
X and Y polarized polariton modes. Therefore, here we 
present the study of only one X polarization. 

Fig|3] (a-f) shows spectrally resolved 2D images 
recorded in momentum space for X-polarized modes, 
each panel corresponding to the energy of each mode 
for the same excitation conditions and spot size as in 
FigOJc). It is seen that the lower energy states are local- 
ized around k = 0. On the contrary, higher energy con- 
densed states have well defined k-vectors different from 
zero. State 1, which is the ground state mode emitting 
at 1.67429 eV, is localized around k = within ±0.6 
/im _1 . The second and third modes (state 2, state 3) 
emitting at energies 1.67445eV and 1.67455eV, respec- 
tively, show increasing spreading in k space. The fourth 
and fifth modes show multiple lobe structure (state 4, 
state 5) and the fifth mode has four well defined sym- 
metrically positioned lobes. The lobes have wave vectors 
(k x ,ky) w (-1, -l)/ira _1 , (k x ,ky) « (1, l)/ira _1 and 
(k x ,ky) « (lj-l^ra -1 ,^,^) « (-1, l)/ira _1 . The 
highest energy mode (state 6) is the weakest and is found 
on a ring in k-space with most of the population being 
at the negative k x side. Simultaneous appearance of con- 
densate density at ±k was observed in numerous other 
k-space images that we recorded (not shown) and also in 
theoretical simulations (see below). The satellites at — k 
may be related to coherent backscattering of the polari- 
tons that condense at +k due to its interaction with the 
local photonic potential. 



FIG. 6: (a)-(f) contrast of the interference pattern of each 
mode. In the blanc regions, the measured intensity is very 
low giving a non denned contrast. 

FigHl shows spectrally resolved images in real space 
recorded for each of the narrow modes for the domi- 
nant X polarization. The lowest energy mode (state 1) is 
strongly localized in real space within an area of ~ 6 x 2 
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FIG. 7: (a),(b) real space correlation pattern. The blanc 
regions are the regions where there is no intensity and thus 
the correlations lead to a high amplitude noise 
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fim 2 , whereas the wavefunctions of the higher energy po- 
lariton states are observed to consist of several maxima 
separated by « 2 — 8 fim and are extended over a length- 
scale 10 — 15 jiva . Moreover all the modes show strong 
spatial overlap, indicating coexistence of the condensed 
phases. 

As it is shown below there is a long range spatial coher- 
ence established for each of the spectrally resolved con- 
densate modes over distances comparable to their size. 
Therefore, the images of the modes in k-space are ex- 
pected to be Fourier transforms of the corresponding 
images in real space [22]. In other words, the k-space 
patterns arise from interference between emission of the 
same energy originating from different regions in real 
space. The fact that the higher energy modes (states 
4-6) have maxima of intensity at k- vectors different from 
zero suggests that there is a large phase difference be- 
tween polar iton fields, which originate from different ar- 
eas (maxima) in the corresponding real space images. 
The k- vectors of the higher energy condensed modes are, 
thus, expected to correspond to a wavelength of the or- 
der of the separation between the maxima in agreement 
to experimental observations. This is further discussed 
in the theory part of this manuscript. 



B. Spatial coherence of spectrally resolved 
condensate modes 

A defining feature of a Bose-Einstein condensate is long 
range spatial coherence. As in previous experiments on 
the same sample, the correlations of the emitted lumi- 
nescence in real space were measured by means of an 
actively stabilized Michelson interferometer in a mirror- 
retroreflector configuration 0, [ll|, [2l| . In order to probe 
the coherence properties of every mode separately, we 
measured the correlations of each individual mode in real 
space. The interferogram at the output of the interferom- 
eter was directed to the entrance slits of the monochro- 
mator and then the real space interferogram was spec- 
trally resolved. The interferogram that we acquired in 
real space allowed us to gain access to the correlations 
for each state. Fig J5] shows the reconstructed interference 
pattern for each of the modes. The extracted contrast of 
each state in real space reveals the coherent content of 
each of the modes. Figj6] shows the contrast for polar- 
ization X, calculated from the interference pattern. The 
existence of highly correlated polariton populations in 
each of the modes, indicates formation of coherence over 
a length comparable to the mode size (5 — 15 fim). The 
lack of correlations for the ground state mode (state 1) 
is not due to the lack of a coherent polariton population 
but rather due to the positioning and strong localization 
of this mode in real space, which results in no overlap 
between the mode and its retroreflected image. 

FigEl shows energy averaged correlations in real space. 
The correlations are observed to be very low and extend 
over only 2 — 3 /am. Such an observation arises from the 



fact that the condensed modes have different momen- 
tum distributions and, thus, the momentum distribution 
of the total polariton emission is quite broad. Therefore, 
spatial interference patterns arising from each of the con- 
densed modes are different and may be expected to cancel 
out in the case of energy- averaged detection. However, 
we note that it is possible to obtain long range spatial 
coherence for energy averaged emission of polariton con- 
densates Q, if one selects a position on the sample for 
which a single polariton condensed mode with a relatively 
narrow momentum distribution dominates above thresh- 
old. In previous experiments [4], attention was always 
paid to choose a position on the sample with the small- 
est possible disorder and long spatial coherence. Here we 
present the results for a generic position on the sample 
with larger amplitude disorder. 

As the final point of this section, we point out that 
multimode structure is a feature that was not revealed 
when the condensate was created on the same sample 
with a multimode Ti:Sapph laser [J], because in that case 
the condensate linewidth was larger than the mode spac- 
ing. The relevant difference between the Ti:Sapph and 
the diode laser is the time scale of intensity fluctuations. 
For the Ti:Sapph laser, the mode spacing is of the or- 
der of 250 MHz, resulting in intensity fluctuations on ns 
timescale. The intensity fluctuations of the diode laser 
are much faster (ps timescale) due to their much larger 
mode spacing 25 GHz). The importance of the inten- 
sity fluctuation time scale can be understood in the fol- 
lowing way. Slow laser intensity fluctuations are followed 
by the exciton reservoir (relaxation time of the order of 
ns), whereas fast fluctuations are averaged out. A fluc- 
tuating reservoir density results in a temporal variation 
of the blue shift and hence a linewidth broadening. This 
explains why the condensate linewidth under Ti:Sapph 
excitation (hundreds of /ieV) is much larger than the one 
under diode laser excitation (tens of fieV) [To| . 



III. THEORETICAL DESCRIPTION 

The simultaneous occurrence of macroscopically oc- 
cupied polariton states at different frequencies is only 
possible out of equilibrium: the condensate frequency 
corresponds to the chemical potential [23], which for an 
equilibrium system is unique and spatially homogeneous. 
Considering the short life time of polaritons (a few ps), 
the non-equilibrium character is not surprising: e.g. the 
thermal velocity does not allow a polariton to travel more 
than a few microns within its life time, a distance that 
is not sufficient to sample the whole excitation area. Po- 
lariton condensates share this non-equilibrium character 
with conventional lasers. 

The crucial difference between a polariton condensate 
and an ordinary laser lies rather on the microscopic than 
on the macroscopic level: a polariton condensate does 
not require population inversion, because the statistics of 
collisions ensures the relaxation of hot excitons, avoiding 
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FIG. 8: (a) Disorder potential landscape used in the theoret- 
ical simulations and (b) energy spectrum from a simulation 
with the mean field model (pQ) and J2}. Multiple condensates 
are found in the simulation. The zero of the energy is the 
same as in (a). Values of parameters used in the simula- 
tions: hg = 0.015 meV/xm 2 , /ry c = 0.5 meV, Hjr = lOmeV, 
hR[n R ] = (meV/im 2 ) x n R hg R = 0, Q = and P/P th = 2. 



the need for population inversion. On the macroscopic 
level on the other hand, there are no fundamental dif- 
ferences and there is a strong analogy between polariton 
condensates in random potentials and lasing in random 
materials [291 ]. 

Given the good spatial coherence of each individual 
condensate mode (see Fig|6]), it is justified to neglect in 
our theoretical description the fluctuations and to study 
the formation of multiple condensates within a mean field 
model [HI, [l6[ . We will adopt here the equations intro- 
duced in Ref. P3] , that were used to explain the effects of 
a finite size pump spot [24| and the appearance of quan- 
tized vortices [11] in polariton condensates. It consists of 



a Gross-Pitaevskii equation including losses and gain for 
the macroscopically occupied polariton field ^( x ) : 



ih 



^Eo — 

+ V ex t 



(r) + %|^(r)| 2 



■Vfe(r)}v(r), (1) 



where E = Huoq and m are respectively the minimum en- 
ergy and the effective mass of the lower polariton branch 
and g > quantifies the strength of repulsive binary in- 
teractions between condensate polaritons. Exciton and 
cavity disorder can be included as an external poten- 
tial term V ex t(r). At the simplest level, the correspond- 
ing gain rate R[ur] can be described by a monotonically 
growing function of the local density ur(t) of reservoir 
excitons in the so-called bottleneck region [25]. At the 
same time, the reservoir produces a mean- field repulsive 
potential Vr(t) that can be approximated by the linear 
expression Vr(t) ~ HgRTiR(r) + HQ P(r), where P(r) is 
the (spatially dependent) pumping rate and qr,G > are 
phenomenological coefficients The GPE equation (pQ) for 
the condensate has then to be coupled to a rate equation 
for n R (r) : 

h R (r) = P(r) - <y R n R (r) - R[n R (r)} |^(r)| 2 , (2) 

where P describes the filling of the reservoir by the non- 
resonant excitation. 

Let us recapitulate the solutions of the equations (pQ) 
and ([2]) in the spatially uniform case V ex t = 0. For low 
excitation density P, a stable steady state is given by 
i/; = and tir = P/jr. The exciton density reaches 
the threshold for polariton condensation when the gain 




FIG. 9: Momentum space distribution of the polariton field 
at the corresponding peak frequencies in Fig. [SJb). 



FIG. 10: Real space distribution of the polariton field at the 
corresponding peak frequencies in Fig. [HJb). 



7 



from stimulated scattering from the exciton reservoir into 
the lower polariton branch equals the polariton loss rate: 
R(nR) — 7c- For pump powers above the threshold P > 
Pth — lR nt R<) the solution without a condensate ip = 
becomes dynamically unstable. The dynamically stable 
solution is then given by the condition that the gain is 
clamped to the losses R[tir\ = 7 C and the condensate 
density is |^| 2 = (P — Pth)/ 7- The condensate wave 
function oscillates at the frequency uj c = g\i(j\ 2 + Vr. 

In the presence of an external potential, this simple 
picture breaks down. In particular, it is no longer guaran- 
teed that only a single frequency appears above the con- 
densation threshold. In the context of polariton conden- 
sation multiple frequency solutions of the related complex 
Ginzburg Landau equation have been found for certain 
types of regular external potentials [261 ]. 

First of all, the photonic potential leads to the local- 
ization of the lower energy polariton modes with small 
k- vectors on a scale of few microns. However, there is 
a continuum of polariton states with higher k-vectors, 
which have kinetic energies above the amplitude of the 
photonic potential disorder. The losses can be overcome 
by the gain for many modes with different energies, which 
may result in simultaneous condensation into both local- 
ized and delocalized states. 

The disorder potential that we used for the theoret- 
ical calculations is shown in Fig. [5(a). The shape of 
the disorder potential is determined from the spatial dis- 
tribution of the low pump intensity photoluminescence. 
The spatial variation of the lowest emission energy fol- 
lows the disorder potential, albeit not perfectly because 
of the zero point kinetic energy that smooths out the 
variations of the potential energy. To compensate for the 
effect of the kinetic energy, we have multiplied the varia- 
tions in the experimentally observed energy landscape by 
a heuristic factor of 3.5 so as to reproduce approximately 
the separation between the condensate frequencies. 

The solution of our model equations (pQ) and ([2]) is 
shown in Figs. [8(b)-[9j As in the experiment, multiple 
frequencies are present in the spectrum. The observation 
of multiple frequencies is consistent with the parameters 
used in the simulations: the fluctuations in the disorder 
potential are larger than the blue shift due to condensate- 
condensate interactions g\ip\ 2 (its maximal value in the 
simulation is 0.3 meV). The coupling between the various 
modes at different energies is therefore too weak to en- 
sure synchronization [26j]. The real and reciprocal space 
distributions at the peaks in Fig. [8(b) are shown in Figs. 
[TOl and [9] respectively. 

First of all, we want to draw the attention to the strong 
qualitative similarities between the theoretical simula- 
tions and the experimental measurements. In both cases, 
the lowest energy states are very localized in real space 
(extended over a few microns). In contrast, the states at 
higher energy are delocalized: they are extended over a 
much larger area, with a diameter of the order of 10 /mi. 
The qualitative correspondence between the experimen- 
tal and theoretical reciprocal space images is also excel- 
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FIG. 11: Real space density of the calculated lowest energy 
microcavity eigenstates that are localized within the excita- 
tion region. Modes A and B and C can be identified with the 
three lowest condensate modes. 



lent. The low energy states form a disk around k = 0, 
while the states at highest energy lie on a ring in momen- 
tum space. Note that all the modes except the lowest one 
do not exhibit inversion symmetry in /c-space. This lack 
of inversion symmetry is also present in the experimental 
momentum distributions above threshold: we have ob- 
served many peaks in the momentum distribution which 
have a satellite at the opposite wave vector. This is likely 
to be related to the phenomenon of backscattering in dis- 
ordered systems. 

Time reversal of the linear Schrodinger equation re- 
quires that the eigenstates in the disorder potential are 
invariant under k — » — k. The condensate modes there- 
fore do not coincide with the linear eigenstates of the 
microcavity. An identification of the condensate states 
with the linear eigenstates (see Fig. [11]) can only be made 
for the lowest condensate state (it corresponds with the 
lowest eigenstate A). The higher energy states, on the 
other hand cannot be identified with a linear eigenstate. 
Instead, the nonlinear dynamics of pumping, dissipa- 
tion and polariton-polariton interactions defines the new 
modes into which the polaritons condense. 

We now discuss what selects the energies and k-vectors 
of the delocalized modes at which polariton condensation 
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state 1 stated 




FIG. 12: Simultaneous representation of the density (hight of 
the surface) and phase (color coding, lighter is higher phase) 
of the polariton condensate modes. Note that the phase dif- 
ference between neighboring maxima is typically large. 

is triggered. To obtain further insight into the nature of 
the condensate states, it is instructive to look simultane- 
ously at the density and phase of the polariton modes. 
The height of the surface of FigHS represents the den- 
sity, where the color represents the phase. Within one 
maximum of the polariton density, the phase is almost 
constant. The phase difference between the peaks is by 
contrast large (of the order of tt). These phase differences 
allow for currents between the different potential minima 
in order to lock to the same frequency [27]. These states 
are reminiscent of the excited ir states that were observed 
with polariton condensates in a periodic instead of a dis- 
ordered potential [28] . Josephson currents are associated 
with the large phase difference between the different con- 
densate regions. They flow from high (light colored) to 
low (dark colored) phase. The large phase difference be- 
tween the polaritons at the different potential minima 
explains why the radius of the momentum distributions 
corresponds to a wavelength of the order of the separa- 
tion between the potential minima. Indeed, as seen from 
the theoretical and the experimental real space images 
of the condensed states (FigJH and FigOU respectively) 
the typical separation d between the intensity maxima 
in real space is about 3 — 4 /am. Therefore, taking into 
account the large phase difference the emission in the 
far-field due to interference is expected to have maxima 
at angles with respect to the normal of the samples 
such as dsm(0) = A/2, where A is the wavelength of con- 
densate emission. These angles correspond to the values 
of k- vectors k = n/d ~ 1/im -1 in agreement with the 
experimental and theoretical k-space images (Fig|5] and 
FiglHl respectively). Note that in the second state (state 
2, Fig fTUlf c)), the flow of polaritons is directed from the 
outside to the center, thereby avoiding in plane losses 



of polaritons out of this condensate mode. This idea of 
macroscopic occupation of states that minimize the losses 
has been put forward in the theory of random lasers [29[ . 
However, in the present system, such a a mechanism does 
not seem to be important: the flow in the third state 
(state 3, FigfTQlc)) goes in the outward direction, yet the 
largest number of polaritons is condensed in this state. 



IV. CONCLUSIONS 

We have presented a detailed experimental and the- 
oretical analysis of a multimode polariton condensate. 
The real and reciprocal space densities of the different 
modes were recorded. We have observed that the lowest 
energy mode is typically localized within a single deep po- 
tential minimum, whereas the higher energy modes are 
delocalized over many potential minima. In reciprocal 
space, the lowest energy mode is centered around k = 0; 
the higher energy modes tend to have a /c-space density 
distributed on a ring with a radius that corresponds to a 
wavelength of the order of the distance between the po- 
tential minima. Due to the breaking of time reversal, the 
condensate modes lack inversion symmetry in reciprocal 
space, although often a peak observed in the momen- 
tum distribution at given k is accompanied by a weaker 
peak at — k, a phenomenon which may be explained by 
the occurrence of coherent backscattering in the disorder 
landscape. Finally we have experimentally demonstrated 
that despite the low long range spatial coherence in the 
multimode state, energy resolved interferometry shows 
good spatial coherence for the individual modes. 

Given the good spatial coherence of the individual 
modes, it is justified to perform the theoretical analysis 
in the mean field approximation. We have used a model 
that describes the polariton dynamics with a generalized 
Gross-Pitaevskii equation coupled to a rate equation for 
the density in the exciton reservoir. With this model, 
we were able to reproduce the main features of the real 
and reciprocal space densities. The analysis of the con- 
densate modes showed that they are in general different 
from the linear eigenstates in the disordered potential, 
but that the lowest energy state is typically very simi- 
lar to the lowest energy eigenstate within the excitation 
spot region. The large phase difference of the polariton 
field between different potential minima explains why the 
radius of the rings in momentum space correspond to a 
wavelength of the order of the separation between the 
potential minima. 
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